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Abstract. Trifluoroacetic acid (TFA; CF 3 COOH) is produced by the degradation of the 
halocarbon replacements HFC-134a, HCFC-124, and HCFC-123. The formation of TFA 
occurs by HFC/HCFC reacting with OH to yield CF 3 COX (X = F or Cl), followed by 
in-cloud hydrolysis of CF 3 COX to form TFA. The TFA formed in the clouds may be 
reevaporated but is finally deposited onto the surface by washout or dry deposition. 
Concern has been expressed about the possible long-term accumulation of TFA in certain 
aquatic environments, pointing to the need to obtain information on the concentrations of 
TFA in rainwater over scales ranging from local to continental. Based on projected 
concentrations for HFC-134a, HCFC-124, and HCFC-123 of 80, 10, and 1 pptv in the year 
2010, mass conservation arguments imply an annually averaged global concentration of 
0.16 jag/L if washout were the only removal mechanism for TFA. We present 3-D 
simulations of the HFC/HCFC precursors of TFA that include the rates of formation and 
deposition of TFA based on assumed future emissions. An established (GISS/Harvard/ 
UCI) but coarse-resolution (8° latitude by 10° longitude) chemical transport model was 
used. The annually averaged rainwater concentration of 0.12 /ig/L (global) was calculated 
for the year 2010, when both washout and dry deposition are included as the loss 
mechanism for TFA from the atmosphere. For some large regions in midnorthern 
latitudes, values are larger, 0.15-0.20 ^tg/L. The highest monthly averaged rainwater 
concentrations of TFA for northern midlatitudes were calculated for the month of July, 
corresponding to 0.3-0.45 /ag/L in parts of North America and Europe. Recent laboratory 
experiments have suggested that a substantial amount of vibrationally excited CF 3 CHFO is 
produced in the degradation of HFC-134a, decreasing the yield of TFA from this 
compound by 60%. This decrease would reduce the calculated amounts of TFA in 
rainwater in the year 2010 by 26%, for the same projected concentrations of precursors. 


1. Introduction 

The long-lived fully halogenated organic compounds are 
being replaced by hydrofluorochlorocarbons (HCFCs) and hy- 
drofluorocarbons (HFCs). Because of removal by reaction with 
OH in the troposphere, smaller tropospheric accumulations 
for these chemicals are expected [World Meteorological Orga- 
nization ( WMO ), 1990; Intergovernmental Panel on Climate 
Change (IPCC), 1990; Rodriguez et al., 1993; Kanakidou et al., 
1995], leading to smaller global ozone depletion potentials and 
global warming indices. Another criterion for the acceptability 
of the substitute HCFCs/HFCs pertains to the impact of these 
compounds or their degradation products on ecosystems. 

Trifluoroacetic (TFA; CF 3 COOH) is produced from the at- 
mospheric degradation of HFC-134a (CF 3 CH 2 F), HCFC-124 
(CF 3 CHFCI), and HCFC-123 (CF 3 CHC1 2 ). It is highly soluble 
in water and expected to be removed from the atmosphere 
mainly by washout. Globally averaged concentrations of TFA 
in rainwater (pg/L) r TFA can be easily estimated from mass- 
balance considerations. Adopting molar yields for TFA of 33% 
yield from HFC-134a (see discussion in section 5.6 for impact 
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of other possible yields) and 100% for HCFC-123 and 124, a 
global rainfall rate of 5 x 10 17 L/yr [Jaeger, 1983] and an 
emission scenario for the HCFC/HFCs provided by D. Hufford 
(1993), the global average r TFA concentrations for the year 
2010 is about 0. 16 pg/L, assuming that washout is the dominant 
loss mechanism for gas phase TFA. 

Local r T FA values deviate from the global averages due to 
inhomogeneities in OH concentrations and precipitation rates, 
as well as the long-range transport of gas phase. TFA and its 
precursors. Evaluation of the impact of these processes re- 
quires models ranging in scales from global to regional. The 
global distribution of TFA in rainwater produced from HFC- 
134a alone has been previously studied with a two-dimensional 
(2-D) model [Rodriguez et al., 1993]. Kanakidou et al. [1995] 
modeled the fate of a series of CFC replacements, including 
TFA precursors, with a 3-D model using monthly average wind 
fields and full chemistry. 

In this paper we present the global distribution of TFA in 
the gas phase and in precipitation produced by the degradation 
of HFC-134a, HCFC-124, and HCFC-123 using a coarse grid 
chemical transport model (CTM), but including synoptic scale 
weather patterns, diurnal boundary layers, and wet and dry 
convective plume mixing processes derived from a general 
circulation model (GCM). Observed climatological values of 
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Figure 1. Degradation pathway for HFC-134a from Franklin 
[1993], Diagram shows mechanisms that lead to the formation 
of TFA. CF 3 radical formed in the end undergoes several other 
reactions, which arc of no importance to the topic discussed 
here. 


cloud cover and rainfall were used to calculate hydrolysis and 
washout rates. On an annual average basis in midnorthern 
latitudes, r TFA concentrations are 0.15-0.20 /ng/L in specific 
grid boxes, but for particular months they are as large as 
0.3-0.45 jug/L. Modeling uncertainties due to assumed liquid 
water and precipitation patterns are also presented. 

No deleterious effects have been identified for the expected 
water concentrations of less than 0.5 jutg/L. Concerns have been 
raised about the long-term accumulation of trifluoroacetic acid 
in certain aquatic environments. Since TFA seems to be re- 
markably stable in the biosphere/hydrosphere, Tromp et al. 
[1995] suggested that accumulation of TFA in certain closed 
aquatic systems over several decades could increase TFA con- 
centrations to factors of 100 times rainwater concentrations. 
The study presented in this paper provides values for rainwater 
concentrations averaged over areas of the order of 800 X 1000 
km. This information can be used to assess the effect of long- 
term accumulation by obtaining more information on the local 
hydrological conditions, atmospheric degradation of TFA pre- 
cursors, and the local deposition fluxes and the local concen- 
trations in rainwater of TFA. Individual aquatic environments 
in a suspected impact area must be studied by using hydrolog- 
ical models with the drainage basin characteristics of the area. 

The degradation pathways for the precursor HFC and 
HCFCs are discussed in section 2. A brief description of the 
chemical model, input data sets used for rainfall and cloud 
cover, and approximations used for calculating cloud and pre- 
cipitation processes are presented in section 3. Section 4 pre- 
sents the results of the simulations. Sensitivity calculations for 
noted uncertainties are discussed in section 5. Conclusions and 
recommendations for further research are found in section 6. 


2. Degradation Pathway for HFC/HCFCs 

Our simulation of TFA concentrations in rainwater requires 
calculation of the net yield of TFA in the degradation of 


HCFC/HFCs, which has been the subject of intense research 
for the past decade. Schemes originally proposed in analogy 
with those of fully saturated hydrocarbons [WMO, 1990] were 
proven correct by numerous kinetic measurements [e.g., 
NASA , 1994]. We review here the degradation of the three 
TFA precursors, HFC- 134a, HCFC-123, and HCFC-I24, 
whose main atmospheric sink is initiated by reaction with OH 
in the troposphere. Stratospheric removal occurs also through 
photolysis and less so by reaction with 0( I D). 

Figure 1 shows our current understanding of the degrada- 
tion pathways for HFC- 134a [Franklin, 1993]. Analogous deg- 
radation mechanisms have been described for HCFC-123 and 
124 [see WMO , 1990]. Consideration of the degradation path- 
ways in Figure I indicates that the TFA rainout is a major 
pathway for final removal of the trifluoroacetyl fluoride, 
CF 3 C(0)F, from the troposphere. The actual yield of TFA will 
depend on the competing photochemical pathways for the deg- 
radation of alkoxy radical (CH 3 CHFO). One option is that the 
complete photochemical scheme be included in a 3-D CTM; 
however, transport will affect the geographical distribution of 
TFA only through its impact on those intermediate products 
whose lifetimes are longer than a few' days in the lower tropo- 
sphere, where most of the rainout occurs. Thus we simplified 
the chemical scheme and were able to carry out a larger num- 
ber of sensitivity studies. 

For the purposes of the TFA simulations the schemes for 
precursor degradation (Figure 1) were simplified to those 
shown in Figure 2, in which only three species with lifetimes 
longer than a few days in the lower troposphere are carried in 
the CTM as prognostic variables. These include the precursor 



Figure 2a. Degradation pathways for HFC- 134a. The reac- 
tion pathway is highly simplified and does not show many 
short-lived intermediate compounds. The three step TFA for- 
mation mechanism hypothesized for the 3-D simulation is 
shown on the left of the figure. Percentages denote the calcu- 
lated molar fraction on each box degraded through the indi- 
cated channel. See text for discussion of a. 
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HCFC/HFC, an intermediate trifluoroacetyl halide of the form 
CF 3 C(0)X (X = F, Cl), and gas phase TFA. Figures 2a, 2b, 
and 2c show the intermediate steps in the formation of TFA 
from HFC-134a, HCFC-123, and HCFC-124. Local photo- 
chemical steady state is assumed for all the other species not 
explicitly included in this scheme. 

2.1. Removal of Precursor HCFC/HFC and Production 
of CF 3 C(0)X 

As indicated, in Figure 2, the removal of the precursors 
occurs through reaction with tropospheric OH. The reaction of 
HFC/HCFC with OH leads to production of CF 3 CXYO 
through multiple and rapid intermediate steps (where X = F, 

Y = H for HFC-134a; X = Y = Cl for HCFC-123; and X = F, 

Y - Cl for HCFC-124). 

P(CF 3 CXYO ,r, t) = * 0 „(r)[OH](r, Of HCFC/HFC] (r, /) 

(i) 

where P(CF 3 CXYO, r, t) (cm -3 s -1 ) denotes the production 
of CF 3 CXYO at each location r and time /; k aH (T) (cm 3 s -1 ) 
is the reaction rate (a function of local temperature (T)) of the 
corresponding HCFC/HFC with OH; and the local concentra- 
tions (cm 3 ) of OH and HCFC/HFC are denoted by brackets. 
The rate constants (fc OH ) adopted for the reaction of HFC- 
134a, HCFC-123, and HCFC-124 with OH are from DeMore et 
al [1992]. 

Particular attention must be paid to the degradation of 
HFC-134a. The radical CF 3 CFHO is produced in the atmo- 
sphere mainly via reaction of CF 3 CHF0 2 with NO. Laboratory 
data reported by Wallington et al. [1996] indicate that the 
CF3CHFO is produced initially in an excited state 



Figure 2b. Same as 2a but for HCFC-124. We assume that 
all CF 3 CClFO is thermalized. 



Figure 2c. Same as 2a but for HCFC-123. We assume that ail 
CF 3 CC1 2 ) is thermalized. 


(CF 3 CHFO*). A fraction [-a undergoes rapid dissociation 
through C-C fission, while the remaining fraction a becomes 
deexcited by thermal collisions. The thermalized alkoxy inter- 
mediate CF 3 CFHO can either form CF 3 COF after losing the 
H atom through reaction with 0 2 or undergo thermal fission of 
the carbon-carbon bond. 

CF3CFHO + 0 2 -» CF3COF + HO, (k 2 ; cm 3 s' 1 ) (2) 

CF3CFHO -> CF 3 + HCOF (fc 3 ; s* 1 ). ‘ (3) 

Taking into account the possibility of prompt dissociation of 
CF3CHFO* in the atmosphere, the molar yield of CF 3 COF 
from degradation of HFC-134a is given by 

k 2 (T)[0 2 ] _ [0 2 ] 

yma a k } (T) + k 2 (T)[0 2 ] a K(T) + [ O. W 

where K(T) (cm -3 ) = (k 3 (T)/k 2 )). The rates k 2 and k 3 have 
been measured by Tuazon and Atkinson [1993] and Wallington 
et al. [1992]. The reported values of Wallington et al. [1992] 
agree with those of Tuazon and Atkinson [1993] when the rate 
expression for K of Wallington et al. [1992] is expressed in a 
pressure dependent format as shown in Table 1. Wallington et 
al. [1996] reported a value of 0.4 for a. We note, however, that 
the calculations presented below assume a value of a = 1.0. 
The impact of the recent value of Wallington et al. is discussed 
in section 5.6. 

By contrast, the alkoxy radicals CF 3 C1 2 0 and CF 3 ClFO, 
produced in the degradation of HCFC-124 (Figure 2b) and 
HCFC-123 (Figure 2c), respectively, undergo chlorine abstrac- 
tion leading to the formation of CF 3 COX [Atkinson, 1989], 
with 
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Table 1. List of Reaction Rate and Hydrolysis Rate Constants Used in the Model 


Species 

Constant 

Reference 

CF,COF 

CF.COF 

CF.COCl 

CFjCOCl 

CF3COF 

H = 3.2 (M/atm) (Henry’s law constant) 
k h = 147 (1/s) (hydrolysis constant) 

H = 2.5 (M/atm) (Henry’s law constant) 
k h = 200 (1/s) (hydrolysis constant) 

YiT \r\ *7(7W hi ok«n. { ».tn«/A\(T)n 2 } 1 

ML A/) , + ‘ 

George et al. [ 1 993] 
George et al. f 1 993] 
George et ai [ 1 993] 
George et ai [1993] 


1 

where Kj,T) - , J8 x ]( ,- 2 5 exp 3 «xw 

,-TA 735 

" 3.28 x 10'°*1 .58 x 10 25 exp(3600/7) 
A/ = air density molecuIes/cnL F c = 0.6 



y i23 - yi24 = i (5) 

We can then express the production of CF 3 C(0)X at each 
model grid point using local steady state 

P(CF 3 C(0)X, r, t) 

= yiicFC/iiFc( r * Of'on(r, OfOHlftr), *) [HCFC/HFC](r, t) (6) 

where HCFC/HFC denotes HFC-134a or HCFC-124 if * = F, 
and HCFC-123 if X = Cl. The yields for HCFC-I23 and 
HCFC-124 is 1 and that for HFC-134a is given in (4). 

2,2. Degradation of CF 3 C(0)X 

Hydrolysis of CF 3 C(0)X in cloud water results in the for- 
mation of TFA and is the main loss process for CF 3 C(0)X. 
Because of the large grid size (approximately 800 km x 1000 
km) in the CTM, the cloud occurrence and hydrolysis are 
necessarily subgrid-scale processes. The rate-limiting step for 
the hydrolysis rate is determined by the slowest of several 
processes, including (1) residence time of air masses inside a 
cloud, (2) diffusion of the gas molecule to the liquid droplet, 
(3) diffusion of the molecule into the liquid, and (4) hydrolysis 
inside the droplet. The analysis of Schwartz [1986] shows that 
the timescales of interest arc the hydrolysis time inside the 
droplets and the residence time of an air parcel inside the 
clouds. 

To address this issue, we adopted the following procedure: 
Lelieveld et ai [1989] have calculated the residence time of an 
air parcel inside a cloud for different cloud types. The resi- 
dence time is defined as the ratio of the depth of the cloud and 
the vertical velocity of air through the cloud for majority of the 
cloud types. The results indicate the residence times range 
between 0.2 and 3 hours for various cloud types. Using this 
definition of residence time, we can approximate the amount 
of air processed by a cloud in a given grid box during a given 
time period A/ as the sum of (1) the fraction of initial amount 
of air within the cloud processed during this time and (2) the 
fraction of additional volume of air which flows through the 
cloud in the time period A t. If the hydrolysis rate is sufficiently 
fast, we can neglect the fact that different air masses will reside 
for different lengths of time within the cloud and obtain the 
following approximate expression for the amount of 
CF 3 C(0)X removed in time step Af: 

A[CFjC(0)X](r, t) = - [CF 3 C(0 )X](r, t) fir, t) 

.(1 - exp ( - K/r, f)A0)(Af/T) (7) 


where Af is the hydrolysis time step used in the model ,/(r, t) 
is the fractional cloud volume in the grid, and K tl is discussed 
below. A hydrolysis time step of 2 hours was used. We adopted 
a value of 122 hours for r as an approximate average of the 
residence times for cumulus and stratus. This r was used for 
the entire model domain, and no differentiation was made 
between upper level, midlevel, and low clouds with respect to 
the residence times. The sensitivity of the calculations to this 
assumed residence time is evaluated in section 5.3. We find this 
formulation to be adequate for the purpose of the present 
calculations and on the basis of our sensitivity tests on the 
timescales of hydrolysis involved. 

The in-cloud hydrolysis first-order rate K d (s 1 ) in (7) is 
calculated as follows: in the absence of mass transfer con- 
straints the fraction of gas phase CF 3 C(0)X that reacts per 
unit time by dissolving and hydrolyzing in the cloud water is 
given by \Kanakidoii et ai, 1993; Freiberg and Schwartz , 1980]. 

K d =fk h HRT (LWC) (8) 

where H is the Henry's constant (A// atm; M is moles dissolved/ 
liter of solution); LWC is the dimensionless liquid water con- 
tent of clouds (grams of water/grams of air); R is the gas 
constant; T is temperature (in °K) and k h is the hydrolysis 
constant in rainwater (s _i ). Solubility and hydrolysis rate con- 
stants for CF 3 COF and CF 3 COCl are given in Table I. The 
factor / denotes a correction factor to compensate for diffu- 
sional resistance inside the droplets, corresponding to values of 
0.5 for CF 3 COCl and 0.65 for CF 3 COF [Kanakidou et al . , 
1995]. The hydrolysis rates calculated by (8) are of the order of 
17min in the presence of clouds. Thus the exponential term in 
(7) tends to zero for the hydrolysis time step considered here 
(2 hours), and the amount of CF3C(0)X becomes directly 
proportional to the cloud cover fraction. 

Cloud fraction and liquid water content of the clouds were 
estimated from the zonally averaged data described by Le- 
lieveld et ai [1989]. This data set contains cloud cover and 
liquid water content of the clouds on a 10° latitude by 100 mbar 
resolution and has been used by Lelieveld and Crutzen [1990] to 
study the effect of liquid phase chemistry on ozone levels in the 
troposphere. The data set used has averaged information for 
the four seasons. The zonally averaged data were mapped on 
to a three-dimensional grid using the climatological monthly 
averaged surface rain data of Shea [1986] as a weighting factor. 
Figure 3 shows the annual average and column-integrated liq- 
uid water. Large-scale features of the liquid water distribution 
are qualitatively consistent with observations [ Prabhakara et 
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al, 1983], but a detailed comparison with other data sets, such 
as ISCCP (International Cloud Climatology Project), is beyond 
the scope of this work. 

We also assumed that the cloud fraction given for a grid box 
is constant throughout the month. Thus there is no cloud free 
period for the entire grid box if the data suggest the existence 
of clouds at the beginning of each month. The usage of surface 
rainfall statistics to map the 2-D fractional cloud coverage and 
liquid water content data onto the 3-D grid seems reasonable, 
since excluding a few locations, persistent cloudiness is indic- 
ative of high monthly average precipitation [Wallace and 
Hobbs , 1977]. There are a few exceptions, where persistent 
nonprecipitating low-level cloud decks are observed, with little 
or no accompanying rainfall. The sensitivity of calculated 
to our assumption of cloud frequency is discussed in section 
5.1. 

We also included photolysis and dry deposition of 
CF 3 C(0)X in the simulations. CF 3 COF is photolyzed slowly in 
the troposphere and has photolysis lifetimes of the order of 
thousands of years [Rattigan et al 1991]. CF 3 COCl produced 
from HCFC-123 (Figure 2c) is photolyzed at a much faster rate 
in the troposphere with photolysis lifetimes of the order of 
weeks [Wallington et al. , 1994]. These photolysis channels do 
not lead to production of TFA. We adopted a dry deposition 
velocity of 0.3 cm/s for CF 3 C(0)X over both land and ocean 
surfaces. 

2.4. Production of Gas Phase TFA in the Atmosphere 

Hydrolysis of CF 3 COX in cloud water produces 
CF 3 C(0)0H (TFA). TFA is a strong acid with a pK of 0.25 
[Exner et al., 1992] and is expected to dissociate immediately 
and form CF 3 C(0)0 . Upon cloud reevaporation this anion 
will probably reevaporale as gas phase CF 3 C(0)OH or remain 
a CF 3 COO~ in residual aerosol. In either case we assume that 
both forms will again dissociate after reincorporation into 
cloud or rainwater. TFA was assumed to undergo full evapo- 
ration as many times as the cloud water itself within each grid 
box and time step. In this case, the production of gas phase 
TFA in the model time interval A t is given by 

P(TFA, r, 0 = - A[CF 3 C(0)X](r, /)/Af (9) 

where A[CF 3 C(0)X] is discussed in the previous section. If 
CF 3 COO“ is incorporated into stable salts which do not fur- 
ther dissociate, the production of gas phase TFA would be 
much smaller than calculated below. 

2.5. Removal of Gas Phase TFA 

TFA is highly soluble (Henry’s constant of > 10 7 M atm 1 
[Wine and Chameides , 1990]) and hence can be expected to be 
rained-out at a rate similar to that of HN0 3 , which has a global 
average washout rate of about 9 days. We parameterized TFA 
washout as a continuous process, using the altitude dependent, 
globally averaged first-order loss rates of nitric acid [Logan et 
al . , 1981]. To account for spatial variations, the specified wash- 
out rates were mapped onto a three-dimensional grid using the 
climatological surface precipitation data of Shea [1986] as a 
weighting factor. This approach is new with this CTM; we 
neglect the coupling between rainout and convective and/or 
large-scale precipitation (see discussion by Balkanski et al 
[1993]). 

Dry deposition of TFA was modeled by assuming a constant 
deposition velocity of 2 cm/s overland and 0.3 cm/s over ocean 
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Figure 3. Annual average column-integrated cloud liquid 
water content (gm/m 2 ). 


surfaces at rates approximately similar to those used for HNQ 3 
[Kotamarthi and Carmichael , 1990]. The deposition velocities 
over the ocean surface are lower as the surface roughness over 
the ocean is smaller, increasing the transfer resistance of the 
gas phase TFA to the ocean surface. Rates for the reactions of 
gas phase CF 3 C(0)OH with OH have been measured [Carr et 
al . , 1994]. Kanakidou et al. [1995] estimated that inclusion of 
this process would decrease the calculated TFA in rainwater 
(r tfa) by at most 5% and much less if a large fraction of the 
reevaporated anion takes the form of a salt. We have not 
included this process in our calculations. 

3. Description of Model and Input Data 

The 3-D CTM developed at GISS (Goddard Institute for 
Space Science)/Harvard/UCI (University of California at Ir- 
vine) was used for simulating the 3-D distribution of the trace 
species of interest. This model uses a 1 year output of wind 
data and convection statistics (every 4 hours) from the GISS 
GCM [Hansen et al., 1983]. This data set is recycled for mul- 
tiple-year runs. The interhemispheric transport in this model 
has been previously tested through comparisons with measure- 
ments of long-lived fluorocarbons, CH 3 CC1 3 and 85 Kr [Prather 
et al . , 1987; Jacob et al . , 1987]. The model reproduced obser- 
vations satisfactorily after corrections were applied to the 
transport scheme to include a horizontal diffusion term to 
improve the north-south hemispheric exchange [Prather et al., 
1987]. Vertical transport by convection was tested through 
comparisons to 222 Rn [Jacob and Prather , 1990]. A coarse-grid 
tropospheric version of the model, with a horizontal resolution 
of 8° latitude by 10° longitude and nine grid points in the 
vertical extending from the surface of 10 mbar was used in the 
present study. There are seven grid points in the troposphere 
and in the stratosphere. 

Simulations of the concentrations of TFA in rainwater were 
carried out by simultaneous calculations of three sets of tracers 
(see Figure 2): (1) the precursor HCFC/HFC, (2) the long- 
lived halocarbonyls of the form CF 3 C(0)X, and (3) gas/ 
particle phase TFA. Production and loss of these tracers are 
incorporated in the model as described in section 2. The pa- 
rameterizations for these processes require input of (1) emis- 
sion scenarios and geographical distribution for precursors 
HCFC/HFCs, (2) global fields of OH for destruction of pre- 
cursors in troposphere, and (3) three-dimensional distributions 
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year 


Figure 4. Emission scenario for the three TEA precursors 
used in the study [Hufford, 1993], 

of liquid water content and precipitation to estimate rates of 
hydrolysis and washout. 

The emissions of precursors between 1985 and 2010 were 
prescribed according to the scenario shown in Figure 4, which 
predicts FTFC-134a to be in wider use than HCFC-124 or 
HCFC-123 during the next decade (D. Hufford, 1993). Other 
emission scenarios suggest that the emissions in Figure 4 could 
easily change by factors of 2 [Kanakidou et al. , 1995]. The 
emissions were distributed around the globe using the electric 
power consumption statistics as a weighting factor following 
Prather et al. [1987]. This emission distribution is similar to the 
chlorofluorocarbon emission inventory of McCulloch et al. 
[1994]. 

A global scale OH field derived for this model, using a 
parameterized chemistry scheme by Spivakovsky et al. [1990], 
was used to calculate the loss of HFC-134a, HCFC-124, and 
HCFC-123 by reaction with OH. The above authors calculated 
this field based on available and assumed global distributions 
of H 2 0, CO, CH 4 , NO,, and 0 3 . The CH 3 CC1 3 lifetime cal- 
culated by using the model is 5.8 years. This lifetime indicates 
that the global OH calculated in this model is about 20% lower 
than the latest recommendations [Prinn et al ., 1995; IPCC , 
1996]. The sensitivity of the calculated r to a 20% change in 
the OH fields is discussed in section 5.5. 

A precipitation climatology based on a historical set of ob- 


servations from the year 1950 to 1979 [Shea, 1986] was used in 
this study to scale the washout rates, cloud cover fractions, and 
liquid water content of the clouds. Data arc available to pro- 
duce the total rainfall contours for each month on a 2.5 lati- 
tude by 2.5 longitude grid. Full longitudinal coverage is derived 
from 45°S to 90 °N. For completeness we appended the zonal 
mean rainfall data from the Jaeger [1976] climatology to the 
Shea data south of 45°S in order to include the effect of 
removal in this region. However, because of uncertainties in 
the rainfall climatology we ignore results south of 60°S and 
north of 80°N. The sensitivity of the calculated midlatitude 
r rFA t° l^e uncertainties in the rainfall amounts from 60°S to 
90°S is discussed in section 5. 

Time dependent calculations were carried out between the 
years 1989 and 2010 separately for HCFC-123, HCFC-124, and 
HFC-134a by varying the surface fluxes of these species ac- 
cording to the emission scenario in Figure 3. The model was 
initialized in 1985 with negligible mixing ratios of the precur- 
sors and degradation products CF 3 C(0)X and TFA. All other 
production and loss processes were kept constant throughout 
the simulation, since predicted changes in overhead ozone and 
tropospheric OH are not expected to change our results sub- 
stantially, and we cannot reliably predict the secondary chem- 
ical effects. 

4. Results 

4,1. Concentrations and Distribution of Precursors 

The calculated global average mixing ratios for HFC-134a, 
HCFC-124, and HCFC-123 between 1985 and 2010 are shown 
in Figure 5. We obtain values of 80 pptv for HFC-134a, 9.8 
pptv for HCFC-124, and 1.1 pptv for HCFC-123 in 2010. Mod- 
el-calculated lifetimes for HFC-134a, HCFC-124, and HCFC- 
123 are 14.4, 6.2, and 1.4 years, respectively. These are within 
the uncertainties of the latest recommendations for these spe- 
cies [IPCC, 1996; Prather and Spivakovsky , 1990], 

The calculated surface distributions of HFC-134a for the 
year 2010 are shown for the month of July in Figure 6. The 
largest gradients are calculated near the equator and close to 
major sources. The southern hemisphere mixing ratios average 
78 ppt during July (southern hemisphere w'inter) and the 
northern hemisphere mixing ratios average about 90 ppt. The 
January mixing ratios in the southern hemisphere (not shown 
here) arc almost uniform at about 73 ppt and above 85 ppt in 
the northern hemisphere. 

Figure 7 shows the calculated surface distribution of HCFC- 
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Figure 5. Global average mixing ratios (pptv) for HFC-134a, 
HCFC-124, and HCFC-123. 



Figure 6. HFC-134a surface distribution for the month of 
July in the year 2010 in units of pptv. Contour levels at 75, 80, 
85, 90, 95, and 110 pptv. 
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Figure 7. HCFC-124 surface distribution for the month of 
July in the year 2010 in units of pptv. Contour levels at 8, 10, 
12, 13, and 15 pptv. 


124 for the year 2010 for the month of July. The southern 
hemisphere mixing ratios are between 10 and 8 pptv for the 
winter months. Figure 8 shows the calculated surface distribu- 
tion of HCFC-123 for the year 2010 for the months of July. 
HCFC-123 has a much shorter lifetime of 1.4 years, producing 
the largest gradients. The southern hemisphere mixing ratios 
are less than 1 ppt for July. 

4.2. Distribution and Fate of CF 3 C(0)X 

Kinetic data [Woffington et al 1992; Taazon and Atkinson, 
1993] indicate that the yields of CF 3 COF from HCFC-124 and 
CF3COCI from HCFC-123 are close to 100% on mole basis, 
corresponding to 87 and 85% on a mass basis for 123 and 124, 
respectively. Figure 9 shows the calculated mass yields of 
CF3COF from HFC-134a as a function of latitude and height. 
Each latitude points represents an average value for all the 
longitudes in that latitude/altitude level. Small yields of (—0.2) 
near the surface increase to larger yields (~0.8) in the upper 
troposphere. The calculated global-average mass yield of 
CF3COF from HFC- 1 34a was 0.37. 

The mass percentages of CF 3 C(0)X lost by the different 
processes are summarized in Table 2. Cloud hydrolysis, the 
channel producing TFA, by far dominates all the CF 3 C(0)X 
degradation mechanisms. Tropospheric lifetimes of 2.6 and 2.1 
days are calculated for CF 3 C(0)F and CF 3 C(0)CI, respec- 
tively. 

Figures 10a and 10b show the surface distribution of 
CF 3 COF and the zonal average latitude-height distribution of 
CF 3 COF, respectively, for the month of July for the year 2010. 



Figure 8. HCFC-123 surface distribution for the month of 
July in the year 2010 in units of pptv. Contour levels at 0.65, 
1.0, 2.0, 3.5, and 4.0 pptv. 



Latitude 

Figure 9. Zonal average mass yield of CF 3 COF from HFC- 
I34a calculated by using the kinetic data of Wattington et al. 
[1992] and a value of a ~ 1 in equation (4). 

The surface distribution of CF 3 COF (Figure 10a) shows higher 
mixing ratios in the tropical regions and near-source locations 
in western Europe and North America. The zonal mean cross 
section shows increasing mixing ratios with altitude. CF 3 COF 
in the lower and middle troposphere is rapidly removed by 
hydrolysis in clouds resulting in lower mixing ratios. Figures 
1 la and 1 lb show the surface distribution of CF 3 COCI and its 
zonal average latitude-height distribution, respectively, for the 
month of July for the year 2010. The surface mixing ratio 
distribution of CF 3 COCI (Figure 11a) again shows features 
similar to those of CF 3 COF, with maximum mixing ratios in 
near-source locations and in the subtropical regions. The zon- 
al-mean cross section shows a similar behavior compared to 
CF3COF. The mixing ratios decrease with altitude in the trop- 
ical and mid-latitude region up to a height of 700 mbar, cor- 
responding to the maxima in cloud cover fraction climatology 
adopted by us. At altitudes from 700 mbar to the tropopause 
levels the mixing ratios increase. At altitudes corresponding to 
the stratosphere in the model, the mixing ratios decrease again, 
as the photolytic loss of CF 3 COCl starts controlling its distri- 
bution. 

4.3. Deposition and Rainwater Concentrations of TFA 

We consider two parameters which could be of interest in 
evaluating the effect of TFA on a given environment. These 
are (1) the deposition fluxes (mass/unit area) of TFA at a given 
location and (2) the concentration of TFA in rainwater (r-y^). 
The deposition fluxes could be a useful input for estimates of 
accumulations of TFA over an extended period of time, par- 
ticularly if we are interested in evaluating the input to a large 
stable water reservoir or soil. Concentrations of TFA in rain- 
water are more pertinent to evaluate the concentrations in 
“seasonal” water bodies produced each year exclusively by 
local rainfall and whose content is wholly determined by yearly 
rainfall and possible year-to-year accumulation [Tramp et aL , 
1995]. 


Table 2. Percent of CF 3 C(0)X (Global Average by Mass) 
Lost Through Different Processes Modeled 


Process 

CF 3 C(0)F 

CF 3 C(0)C1 

Photolysis 

0 . 1 % 

2 . 6 % 

Surface deposition 

2.9% 

3.5% 

Hydrolysis 

97% 

93.9% 
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Figure 10a. Surface distribution of CF 3 COF (1.0E-3 pptv) 
for the month of June in the year 2010. Contour are at 5, 10, 
25, 50, 75, 90, 100, and 110. 
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Figure 11a. Surface distribution of CF 3 COCl (l.OE-3 pptv) 
for the month of July in the year 2010. Contours ar at 1, 3, 5, 
10, 20, 30, and 40. 



Fluxes of TFA deposited by the washout process, accumu- 
lated over the period of simulation (1989-2010) in mg/m 2 are 
shown in Figure 12. Since most of the CF 3 COX is formed from 
HFC/HCFCs in the tropical regions and CF 3 COF is converted 
to TFA with a very short lifetime (local lifetimes of the order 
of hours), production of TFA maximizes in the tropics. Typical 
washout time constants in the tropics are 3 days. These time 
constants are much smaller than typical meridional transport 
time constants (weeks to months), and most of the TFA pro- 
duced in the tropical regions is deposited in these regions. 
Figure 13 shows the cumulative amount of TFA deposited due 
to both the wet and the dr^ deposition processes, from 1985 to 
2010. Comparison of Figure 13 and 12 indicates that our cal- 
culated dry deposition fluxes are of the order of 20% of the wet 
deposition fluxes. 

Figure 14 shows the calculated annual average r I FA concen- 
trations (pg/L) for the year 2010. The rainwater concentration 
of TFA is defined as the ratio of total TFA deposition by 
washout over a grid box (liters/month) to the rainfall received 
in the grid box (liters/month). Thus r XFA concentration may be 
inversely related to the rainfall received in a particular grid 
box. The maximum r TFA concentrations arc indeed found in 
the region close to maximum TFA production and receiving 
the least amount of rainfall. Region in southwestern United 
States and southern Europe have concentrations of up to 0.2 
pg/L. However, as discussed below, the obtained distribution is 
not simply given by a balance of local production and washout 
of TFA; transport of precursors and gas phase. TFA also plays 
an important role. The global average r TFA is 0.12 pg/L , and 


the subtropical and regions of western Asia have the highest 
r TFA concentrations of up to 0.3 pg/L. The highest concentra- 
tions at any given location in the northern hemisphere were 
obtained for the month of July 2010: 0.45 pg/L at locations in 
southern Europe, 0.3 pg/L for most of North America and 
Europe, and up to 0.4 pg/L for parts of North America and 
southern Europe. 

Because of the fairly long meridional transport time con- 
stants, 2-D model calculations have obtained concentrations of 
TFA in rainwater which can be approximately obtained by 
local balance between production and loss within a latitude 
band [Rodriguez et al , 1993]. This is not the case for the zonal 
asymmetries that show up in this CTM calculation. The annual 
average ratio of maximum r TFA concentration obtained in a 
grid box to the global average r XFA concentration of 0.12 pg/L 
for the year 2010 can be obtained from Figure 14. The maxi- 
mum asymmetry is of the order of 2.5 and is obtained in 
western Asia. Concentrations in most of North America and 
southern Europe are up to a factor of 1.5 higher than the 
global average. In southern United States the r TFA concentra- 
tions are up to a factor 2 higher than the global average 
concentrations. As a reference, Figure 15 shows the annual 
average rainfall amounts used in the simulation, scaled by the 
global mean and then inverted. In the absence of transport we 
would expect that the asymmetries in r TFA would correlate to 
those in Figure 15. This is illustrated in Figure 16, which shows 
the calculated r TFA concentrations asymmetries for a case in 
which the transport of TFA and CF 3 COX was switched off. 
This results in r T -p A asymmetries which arc in much better 



Figure 10b. CF 3 COF mixing ratios (l.OE-3 pptv) averaged 
over a latitude band for the month of July in the year 2010. 



Figure lib. CF 3 COCl mixing ratio (1.0E-3 pptv) averaged 
over a latitude band for the month of July in the year 2010. 
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Figure 12. Cumulative amount of TFA deposited by wet 
processes by the year 2010 at the surface in mg/m 2 . 

agreement with the rainfall asymmetries. Thus inclusion of 
transport substantially “smoothes out” the asymmetries in r TFA 
expected from simple consideration of asymmetries in rainout 
and OH fields. 

The competition between transport and local production 
removal in determining r TFA has an even larger impact on the 
expected variability of r TFA in a given geographical region. The 
r T1A concentrations calculated for the year 2010 in all the grid 
boxes of the model and for all the months of the year (36 
longitude x 24 latitude X 12 months) are presented in the 
form of a box plot in Figure 17. The data were binned by the 
amount of rainfall received in mm/month. The vertical axis 
shows the concentrations of r ThA in ju,g/L. Although the me- 
dian r, FA is fairly constant around 0.12 /xg/L for all the bins, 
the range in the r TFA concentrations at low levels of rainfall is 
larger compared to the bins with larger rainfall amounts. This 
reflects primarily the impact of transport when local removal 
time constants become small. We can thus expect the largest 
uncertainty in the calculation in regions, where many of the 
seasonal water bodies of concern are located. In these cases, 
resolution of synoptic variability in the transport is necessary 
for a better estimate. 

5. Discussion of Uncertainties 

Our approximations for both in-cloud hydrolysis and wash- 
out are very simplistic, assuming an average, continuous rate 


0O*H| 
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Figure 13. Cumulative amount of TFA deposited by wet and 
dry processes by the year 2010 at the surface in mg'm 2 . 
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Figure 14. Annual average TFA in rainwater (/xg/L) for the 
year 2010. 

over the grid box, where in actuality these arc sporadic events. 
A better approach would need to consider, for example, cou- 
pling to moist convection events and large-scale rain [Balkan- 
ski et al . , 1993]. Thus our calculations disregard the coupling 
between clouds, rainout, and transport. We have carried out 
limited simulations to study the sensitivity of our results to our 
assumption of continuous hydrolysis and washout rates versus 
episodic cloud and rain events as w T ell as in-cloud residence 
time. By varying the cloud cover, residence time, and washout 
we can gain information as to the magnitude of the error 
introduced in this decoupling. Finally, we also consider the 
sensitivity to the large uncertainty in the removal rates at high 
latitudes as well as to the adopted concentrations of OH. 

5.1. Sensitivity of TFA Production to Intermittent 
Cloud Cover 

Observations indicate that on a global average, a given lo- 
cation on land can go completely cloud free between 17% and 
19% of the time. Over oceans the completely cloud free time 
is observed to drop to 3% of the total time [ Warren et al ., 1986, 
1988]. We performed numerical experiments, with cloud cover 
switched on for 24 hours every 4 days (75% cloud free) and 
every 7 days (85% cloud free). Loss frequencies for CF 3 C(0)X 
were, however, adjusted to yield the same average frequency 
over the whole period of time. All the simulations discussed 
below were performed for CF 3 COF produced from HFC-134a 
alone and for the year 2010 only. 



0 .5 1.0 2.0 5.0 10.0 99.0 


Figure 15. Annual average asymmetries in rainfall amounts. 
Figure show ratios of global average/grid average amounts. 
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Figure 16. Annual average r XFA asymmetries for the year 
2010 from global average concentrations for the “no transport” 
case. 


In the subtropical regions, TFA production increases by an 
average of 30% with local increases of up to 50%. TFA pro- 
duction at middle and high latitudes changed by only a few 
percent. Thus the effect of intermittent cloud cover on annual 
average TFA production from the hydrolysis of CF 3 COF can 
be expected to result in a maximum of 30% increase on a 
regional scale. 

5.2. Sensitivity of TFA Washout to Intermittent Rain 

The effect of intermittent rainfall on washout rates was 
investigated in detail by Giorgi and Chameides [1985], Thomp- 
son and Cicerone [1982], and Stewart et al [1983] employing 
one-dimensional models. In general, it was noted that the 
increase in frequency of the rainfall shortens the washout life- 
time of the species, and for the case of a highly soluble species, 
the washout rate is proportional to the rainfall frequency. Here 
rainfall frequency refers to the number of rainfall events oc- 
curring in a given time. The higher the frequency, the larger 
the number of rainfall events that occur, eventually leading to 
th e con tinuous rainfall and continuous washout. 

We studied the ratio of r XFA obtained for the continuous 
rainfall case to a case where rain occurs once every 7 days for 
one full day. This is an extreme case. The washout rate and rain 
amount for the intermittent rainfall cases were adjusted such 
that the total rate for a particular month remains the same for 
all the cases. In locations (Indian Ocean and eastern Pacific 
regions) with the highest calculated r XFA , the new r TFA goes up 
by as much as 30%. Increases of up to 20% are obtained in the 
western United States and up to 25% in the southwestern 
United States and southern Africa and China. 

5.3. Sensitivity of r, , A to in-Cloud Residence Time 

We also examined the sensitivity of our calculations by in- 
creasing the in-efoud residence time (with a value of 1.2 hours) 
to 3 hours and decreasing it to 0.8 hours. The residence times 
will have 10-20% effect on r XFA concentrations in locations 
receiving low rainfall (subtropical Asia, Africa, southwestern 
United States, and the northwestern location of South Amer- 
ica). The effect on regions receiving high amounts of rainfall is 
less than 15%, and changes of the order of 1-5% were ob- 
tained over the rest of the domain. 


5.4. Sensitivity of r, , A to Washout South of 60°S 

The washout rales from 45°S to 90 C S were calculated on the 
basis of the climatology of Jaeger [ 1976]. However, values south 
of 60°S are very uncertain, in particular since it is difficult to 
distinguish between rain and snow. The influence of increased 
washout (rainfall) in this zonal region on r XFA concentrations 
were examined by increasing the washout rate beyond 60°S by 
a factor of 2. Increasing rainfall in the southern latitudes will 
decrease r XFA by 10% in some locations to the south of 60°S 
but has no significant effect over the rest of the model domain. 

5.5. Sensitivity of r TFA to OH Concentration 

We also examined the sensitivity of our calculations to a 
20% increase in the OH fields used for calculating the degra- 
dation of the TFA precursor HFC-134a. Results indicate that 
the lifetime of HFC-134a decreases by 16%, and the global 
average HFC-134a mixing ratios is reduced by 6% from the 
base case discussed in section 4. The global average r XFA in- 
creases by 12%. Because of an increase in dry deposition over 
arid subtropical regions of Africa and Asia compared to the 
base case, local increases of the order of 2-10% are obtained 
in these regions. 

5.6. Impact of Reduction in the TFA Yield From HFC-134a 

As discussed above, recent experimental results reported by 
Wattington et al. [1996] suggest that the reaction of CF 3 CFH0 2 
with NO yields primarily vibrationally excited CF 3 CHFO rad- 
icals. These radicals will undergo prompt thermal fission, un- 
less they are thermalized by collisions with ambient molecules. 
The experiments suggest that about 40%> of the excited 
CF 3 CFH0 2 is thermalized. Since reaction with NO is expected 
to be primarily the channel for degradation of CF 3 CFH0 2 , 
these results suggest that the mass yield of CF 3 C(0)F from 
HFC-134a is reduced from 37 to 15%. 

The impact of these reductions on the total amount of TFA 
deposited will depend on the specific scenario adopted for the 
relative emissions of HCFC-123 and 124 and HFC-134a. For 
the scenario adopted in these calculations we estimate that the 


< 

0.8 

0.6 

o 


3 


t 

0.4 

zL 

0.2 



mm of rain 


Figure 17. Boxplot of r XFA versus rainfall amounts for the 
year 2010. Data selected for the plot are from the last year of 
the simulation and all the grid boxes. There is one data point 
per month per grid box. The bar in the middle of the box is the 
median of the data and the boxes represent the middle 50% of 
the data binned in the grid. The whiskers represent the accept- 
able range of the remaining data beyond that located inside the 
box. Outliers are denoted by points. 
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calculated-global-averaged TFA concentrations in rainwater 
will decrease by 26%, i.e., from 0.12 to 0.09 jug/L. 

6. Conclusions 

We have carried out simulations of the degradation of HFC- 
134a and HCFC-123 and HCFC-124, with emphasis on the 
production of TFA and its concentration in rainwater. The 
complicated degradation schemes for these precursors have 
been considerably simplified, taking advantage of the fact of 
the short chemical time constants of most intermediate species. 
We tested the sensitivity of the calculated rainwater concen- 
tration of TFA (r XFA ) to various parameterizations and ap- 
proximations adopted in the model. Our limited sensitivity 
studies suggest that these approximations would not have a 
greater than 30% effect on the calculated Ttfa- 

Annual average concentrations in the range of 0.15-0.2 
jig/L were obtained for the year 2010 in the northern midlati- 
tudes. The calculated r TFA concentrations in July 2010 in parts 
of North America and southern Europe were found to be 
0.3-0.45 {JigfL. Inclusion of a reduced yield for TFA formation 
from HFC-134a, suggested by Wallington et al. [1996], will 
reduce these numbers 26%. Dry deposition fluxes are on av- 
erage about 20% of the wet deposition and about equal for and 
regions. In the absence of soil sinks, dry deposition could also 
contribute to accumulations in groundwater and seasonal wa- 
ter bodies. Assessment of this contribution would require a 
detailed hydrological model for the catch basins of interest. 

Our simulations show the crucial role played by long-range 
zonal transport of gas phase TFA and CF 3 C(0)X in smoothing 
out the chemical impact of asymmetries in the OH and rainfall 
fields. Transport effects become important when advective 
time constants become comparable to those of other transfor- 
mation and deposition processes, a condition that is true most 
everywhere for gas phase TFA. Thus in spite of the short time 
constants for hydrolysis and washout, simulation of the TFA 
distribution does require use of three-dimensional models. 
Transport processes are particularly important in regions, 
where removal rates are slower. Our results show that these 
regions exhibit the largest variance in the calculated Arp/v 

The results presented here can be tested as measurements of 
TFA, and its precursors become available. The calculated mix- 
ing ratios of HCFC-124 and HCFC-123 for NH are 1.96 and 
0.21, and global averages are 1.56 and 0.14 pptv, respectively, 
for July 1995. For the 1995 conditions the contributions of the 
three precursors to the total r TFA in western Europe (0.014 
/ag/L) break down as follows: 57% from HCFC-124, 27% from 
HCFC-123, and 16% from HFC-134a. Thus in the short term 
the greatest impact on could be from shorter-lived 

HCFC-124 and HCFC-123, and it is advisable to monitor their 
emissions. 

Measurements of HFC-134a, starting from the year 1990 to 
1995, were reported by Montzka et al. [1996] and Oram et al. 
[1996]. Montzka et al. [1996] report a mixing ratio of 2.33 pptv 
in the NH, 1.13 in the SH, and a global average mixing ratio of 
1.73 pptv for the beginning of July 1995. The model-calculated 
mixing ratios for this period are 3.4 pptv for the NH, 1.48 pptv 
for the SH, and a global average of 2.44 pptv. This discrepancy 
is a reflection of the fact that this study used projected rather 
than actual emissions for this compound for the period of 
interest. 

Measurements of TFA in natural water bodies and rainwa- 
ter in the western United States were reported by Zehavi and 


Sieber [1996]. It is difficult to compare their measurements for 
natural waters, since we have not carried out a detailed study 
of the local chemistry, transport, and hydrology. Measured 
rainwater concentrations range from 0.03 to 0.09 pgTL, a factor 
of 2-6 larger than our calculated values for the western United 
States in 1995 (0.015 jug/L). Measurements of r TFA have been 
made at two locations in western Europe [Frank et al , 1996]. 
The measured r XFA is in the range of 0.025-0.28 jag/L for the 
year 1995 at Bayreuth, Germany, and in the range of 0.04-0.08 
jug/L at Zurich and Alpthal, Switzerland. The mean r TFA for all 
the measurements is 0.1 pg/L. The measured gas phase con- 
centrations of TFA are of the order of 0.01 pptv. The gas phase 
TFA mixing ratios calculated in the present study are in the 
range of 0.007-0.005 pptv for western Europe for the year 
1995. The calculated annual average r XFA for the year 1995 is 
0.014 pg/L in western Europe. These numbers would have to 
be rescaled when measurements of the precursors are avail- 
able. In the case of HC 134a, this would imply a contribution 
of 0.0015 pg/Liter, based on the measurements of Montzka et 
al. [1995]. Although no measurements for HCFC 123 and HFC 
124 are available, it is doubtful that the actual concentrations 
are higher than the calculated concentrations. 

The higher r XFA concentrations measured could be due to 
additional unknown sources as suggested by Frank et al. [1996], 
We note, however, that given the limited number of samples 
collected and the large variability of the measurements, the 
reported concentrations cannot be considered a true mean for 
the region (i.e., our coarse grid box). Great caution should be 
exercised in comparing the calculated grid average concentra- 
tions and global average concentrations of r XFA to a small 
number of local measurements. Given the sensitivity of the 
calculated present-day r XFA to the amount HCFC-124 and 
HCFC-123 emitted, it would be of great value to measure these 
compounds in the atmosphere. A monitoring program for the 
concentrations of TFA and its precursors could help ascertain 
the actual variability of the rainwater concentrations. 
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